The Behaviour of Varying- Alpha Cosmologies 
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We determine the behaviour of a time-varying fine structure 'constant' a(t) during the early and 
late phases of universes dominated by the kinetic energy of changing a(t), radiation, dust, curvature, 
and lambda, respectively. We show that after leaving an initial vacuum-dominated phase during 
which a increases, a remains constant in universes like our own during the radiation era, and then 
increases slowly, proportional to a logarithm of cosmic time, during the dust era. If the universe 
becomes dominated by negative curvature or a positive cosmological constant then a tends rapidly 
to a constant value. The effect of an early period of de Sitter or power-law inflation is to drive a to 
a constant value. Various cosmological consequences of these results are discussed with reference to 
recent observational studies of the value of a from quasar absorption spectra and to the existence 
of life in expanding universes. 



I. INTRODUCTION 

One of the problems that cosmologists have faced in 
their attempts to assess the astronomical consequences 
of a time variation in the fine structure constants, a, has 
been the absence of an exact theory describing cosmolog- 
ical models in the presence of varying a. Until recently, it 
has not been possible to analyse the behaviour of varying- 
a cosmologies in the same self-consistent way that one 
can explore universes with varying G using the Brans- 
Dicke or more general scalar-tensor theories of gravity. 
However, we have recently extended the generalisation 
of Maxwell's equations developed by Bekenstein so that 
this can be done self-consistently. In a recent letter |l]] 
we reported numerical studies of the cosmological evolu- 
tion of varying-a cosmologies with zero curvature, non- 
zero cosmological constant, and matter density match- 
ing observations. They reveal important properties of 
varying-a cosmologies that are shared by other theo- 
ries in which 'constants' vary via the propagation of a 
causal scalar field obeying 2nd-order differential equa- 
tions. Their structure can be compared with that of 
varying speed of light (VSL) theories developed in refs. 

and with Kaluza-Klein like theories in which con- 
stants like a vary at the same rate as the mean size of 
any extra dimensions of space |2q| . 

Recent observations motivate the formulation and 
detailed investigation of varying-a cosmological theo- 
ries. The new observational many-multiplet technique 
of Webb et al, ||, pQ] , exploits the extra sensitivity 
gained by studying relativistic transitions to different 
ground states using absorption lines in quasar (QSO) 
spectra at medium redshift. It has provided the first 
evidence that the fine structure constant might change 
with cosmological time iPHTll. The trend of these re- 
sults is that the value of a was lower in the past, with 



Aa/a = -0.72 ± 0.18 x 10~ 5 for z « 0.5 - 3.5. Other in- 
vestigations have claimed preferred non-zero values of Aa 
< to best fit the cosmic microwave background (CMB) 
and Big Bang Nucleosynthesis (BBN) data at z « 10 3 



and z 10 10 respectively [|L2||13|, but these need to be 
much larger than those needed to reconcile the observa- 
tions of 

In this paper we present a detailed analytic and numer- 
ical study of the behaviour of the cosmological solutions 
of the varying theory presented in H. We shall confine 
our attention to universes containing dust and radiation 
but analyse the effects of negative spatial curvature and 
a positive cosmological constant. Extensions to general 
perfect-fluid cosmologies can easily be made if required. 



II. A SIMPLE VARYING- ALPHA THEORY 

The idea that the charge on the electron, or the fine 
structure constant, might vary in cosmological time was 
proposed in 1948 by Teller, |14|, who suggested that a oc 
(lni) -1 was implied by Dirac's proposal that G a t" 1 
and the numerical coincidence that a" 1 ~ ln(hc/Gmp r ), 
where m pr is the proton mass. Later, in 1967, Gamow 
|l5| suggested a oc t as an alternative to Dirac's time- 
variation of the gravitation constant, G, as a solution 
of the large numbers coincidences problem but in 1963 
Stanyukovich had also considered varying a, jl(|, in this 
context. It had the advantage of not producing a terres- 
trial surface temperature above 100 degrees centigrade 
in the pre-Cambrian era when life was known to exist. 
However, this power-law variation in the recent geologi- 
cal past was soon ruled out by other evidence. 

There are a number of possible theories allowing for 
the variation of the fine structure constant, a. In the 
simplest cases one takes c and Ti to be constants and at- 
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tributes variations in a to changes in e or the permittivity 
of free space (see || for a discussion of the meaning of this 
choice). This is done by letting e take on the value of a 
real scalar field which varies in space and time (for more 
complicated cases, resorting to complex fields undergo- 
ing spontaneous symmetry breaking, see the case of fast 
tracks discussed in |jj). Thus eo — > e = eoe(a; M ), where e 
is a dimensionless scalar field and eo is a constant denot- 
ing the present value of e. This operation implies that 
some well established assumptions, like charge conserva- 
tion, must give way ^3|. Nevertheless, the principles of 
local gauge invariance and causality are maintained, as is 
the scale invariance of the e field (under a suitable choice 
of dynamics). In addition there is no conflict with local 
Lorentz invariance or covariance. 

With this set up in mind, the dynamics of our theory 
is then constructed as follows. Since e is the electromag- 
netic coupling, the e field couples to the gauge field as eA^ 
in the Lagrangian and the gauge transformation which 
leaves the action invariant is eA^ — > eA^ + x.fi, rather 
than the usual A^ — + A^ +x,/x- The gauge- invariant elec- 
tromagnetic field tensor is therefore 

f»u = \ (M„), M - M„),„) , (i) 

which reduces to the usual form when e is constant. The 
electromagnetic part of the action is still 

S em = - I '^x^f—gF^F^. (2) 

and the dynamics of the e field are controlled by the 
kinetic term 

S. = -~J*xV=9^r, (3) 

as in dilaton theories. Here, I is the characteristic length 
scale of the theory, introduced for dimensional reasons. 
This constant length scale gives the scale down to which 
the electric field around a point charge is accurately 
Coulombic. The corresponding energy scale, hc/l, has 
to lie between a few tens of MeV and Planck scale, 
~ 10 19 GeV to avoid conflict with experiment. 

Our generalisation of the scalar theory proposed by 
Bekenstein described in ref. Q includes the gravita- 
tional effects of ip and gives the field equations: 

= 8nG (T™ atter + T* v + T^e~^) . (4) 

The stress tensor of the ip field is derived from the la- 
grangian — — ^d^ipd^ip and the ip field obeys the 
equation of motion 

□V> = --e-^Lem (5) 

UJ 

where we have defined the coupling constant u> — (he) /I 2 . 
This constant is of order ~ 1 if, as in fit] , the energy scale 



is similar to Planck scale. In order to make quantita- 
tive predictions we need to know how much of the non- 
relativistic matter contributes to the RHS of Eqn. (||). 
This is parametrised by ( = p em / Pmatter- In @, C was 
estimated to be around 1%. However, if we choose to 
model the proton as a charged shell of radius equal to 
the estimated proton radius, the fraction would be lower, 
approximately 0.19% . Also, the value of £ needs to be 
weighted by the fraction of matter that is non-baryonic, 
a point ignored in the literature |lj],|l8|]. Hence, the to- 
tal ( depends strongly on the nature of the dark matter. 
BBN predicts an approximate value for the baryon den- 
sity of Q B m 0.0125/ig 2 , or Q B » 0.03% with a Hubble 
parameter of ho « 0.6. Since we believe the total mat- 
ter density to be Q m w 0.3, this would mean that only 
about 1/10 of matter is baryonic and couples to changes 
in e. Thus, we should assume values for £ ranging from 
0.02% to 0.1%, if the cold dark matter is allowed to have 
little or no electrostatic Coulomb component. If this is 
not true, then £ could have a much higher value. 

We should not confuse this theory with other similar 
variations. Bekenstein's theory fjlj does not take into 
account the stress energy tensor of the dielectric field in 
Einstein's equations, and their application to cosmology. 
Dilaton theories predict a global coupling between the 
scalar and all other matter fields. As a result they pre- 
dict variations in other constants of nature, and also a 
different dynamics to all the matter coupled to electro- 
magnetism. An interesting application of our approach 
has also recently been made to braneworld cosmology in 




III. THE COSMOLOGICAL EQUATIONS 

Assuming a homogeneous and isotropic Friedmann 
metric with expansion scale factor a(t) and curvature 
parameter k in eqn. (Q), we obtain the field equations 
(c=l) 

/ a\ 2 8nG f , > . „ ,,, r „ ,. u) ; 9 \ 

(a) = T" \ Pm ( + C ° XP [ + Pr ° XP [ + 2^ ) 

where A is the cosmological constant. For the scalar field 
we have the propagation equation, 

^ + 3HiP= -exp[-2^]CPm, (7) 

where H = a/a is the Hubble expansion rate. The con- 
servation equations for the non-interacting radiation, and 
matter densities are 

p' m + 3Hp m = (8) 

p r + AHp r = 2ljjp r . (9) 
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and so p m oc a~ 3 and p r e~ 2 ^ oc a -4 , respectively. If ad- 
ditional non-interacting perfect fluids satisfying equation 
of state p= (7 — l)p are added to the universe then they 
contribute density terms p oc a -37 to the RHS of eq.(^) as 
usual. This theory enables the cosmological consequences 
of varying e, to be analysed self-consistently rather than 
by changing the constant value of e in the standard the- 
ory to another constant value, as in the original proposals 
made in response to the large numbers coincidences (see 
ref. |2(J for a full discussion). 

We have been unable to solve these equations in gen- 
eral except for a few special cases. However, as with the 
Friedmann equation of general relativity, it is possible to 
determine the overall pattern of cosmological evolution 
in the presence of matter, radiation, curvature, and pos- 
itive cosmological constant by matched approximations. 
We shall consider the form of the solutions to these equa- 
tions when the universe is successively dominated by the 
kinetic energy of the scalar field ip, pressure- free matter, 
radiation, negative spatial curvature, and positive cos- 
mological constant. Our analytic expressions are checked 
by numerical solutions of (||) and ([?]). 



A. The Dust-dominated era 

We consider first the behaviour of dust-filled universes 
far from the initial singularity. We assume that k = = 
A = py, so the Friedmann equation (|^) reduces to 



8ttG 



( Pm (l + Cexph2^]) + ^ 2 ) , (10) 



and seek a self-consistent approximate solution in which 
the scale factor behaves as 



a — t 



2/3 



dt 



where 



(ipa 3 ) = Nexp[~2^} 



2C 

N = ^p m a 3 

UJ 



(11) 
(12) 

(13) 



is a positive constant. If we put 
x = \n(t) 

then ([[2]) becomes 

?//' + tp' = 7Vexp[-2?/>] 



(14) 



with N > and ' = d/dx. This equation has awkward 
behaviour. For any power-law behaviour of the scale fac- 
tor other than (|llf ) a simple exact solution of (|lj) exists. 
However, the late-time dust solutions are exceptional, re- 
flecting the coupling of the charged matter to the varia- 
tions in ip, and are approximated by the following asymp- 
totic series: 



ip = - \n[2Nx] + ^2 a n x~ n 



(15) 



n=l 



To see this, substitute this in the evolution eqn. (|14| ) for 
ip then it becomes: 



n— 1 



^exp[-2^a n2 ;-"] (16) 



Now we can pick the a n to cancel out all the terms in 
x~ r , r > 2 on the left-hand side. This requires 

02 = 0,1 = — 2' ° 3 = 2 ° 2 ' 0,4 = = 3 x 2a2, etc 



hence 

V- -„ l r l 1 2 

=-2 { x + ^ + ^ 

n—1 

2x3 2x3x4 



(r-1)! 



2T x° 
all that is left of the eqn. (h&j is 



1 1 

2x ~ 2x 



exp[-2 Y] a n x n ] -> 7T 
^— ' 2x 



as x — ► 00. So, at late times, as x = ln(i) becomes large, 
we have 

2x3 2x3x4 (r-1)! 

also, since a = exp[2tp] we have, as t — > 00 

1 1 2 



a = 2N]n(t) x exp[ 



ln(t) (ln(t)) 2 (ln(t)) 3 
2x3 2x3x4 (r-1)! 



(ln(t)) 4 (ln(t))6 
So, to leading order, we have 



a ~2iVm(*)exp[-— ] 



(18) 



(19) 



The non-analytic exp[l/a;] behaviour shows why the 
eqn. (^t[), despite looking simple, has awkward be- 
haviour. We can simplify the asymptotic series ( [18] ) a bit 
further because we know from the definition of the loga- 
rithmic integral function li(x) — J Q dt/hi(t) = Ei[ln(x)], 
that as x — > 00 



li{x) ~ exp[x] J2 -^T 



(20) 



n=0 ' 
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so the series we have in ( |l7| ) in {..} brackets is 

oo / _ j\| 

r ~ ex P [-x]U(exp[x]) (21) 

r— 1 

and so asymptotically 

^ = \ ln[2Nx] - i exp [-x]ii(exp[x]). (22) 
Hence, as t — ► oo, 

</> = ~ ln[2JVln(*)] - i K(t) = i ln[2iVln(f)] - ~Eipn(*)] 



2t 



and so asymptotically, 

a = exp[2tp] = 2N exp[-t~Hi(t)]hit. 



(23) 



(24) 



This asymptotic behaviour is confirmed by solving equa- 
tions (g-||) numerically for p m 3> p r , p^. By using a range 
of initial values for ip we produce the plot in fig (|l|), in 
which the asymptotic solution is clearly approached. 

We need to check that the original assumption of a = 
t 2 / 3 in the Fricdmann eqn. (|J) is self consistent. The 
relevant terms are 



UJ 



p m (l + Cexp [-2^]) + |V 2 



(25) 



The exp[— 2ip] — a 1 falls off as t — > oo so the 
p m (1 + Cexp [—2^]) oc a~ 3 term dominates as expected. 
For the kinetic term ip we have 



(26) 



and so again the ip 2 term falls off faster than t~ 2 as t — > oo 
and the a = t 2 / 3 behaviour is an ever-improving approx- 
imation at late times. If we examine the form of the 



solution (24) we see that a always increases with time 
as a logarithmic power until it grows sufficiently for the 
exponential term on the right-hand side of (Q) to affect 
the solution significantly and slow the rate of increase by 
the series terms. The rate at which a grows is controlled 
by the total density of matter in the model, which is di- 
rectly proportional to the constant N, defined by eqn. 
(|l3|). The higher the density of matter (and hence N) 
the faster the growth in a. However, because of the log- 
arithmic time- variation, the dependence on p m ,u), and £ 
is weak. The self-consistency of the usual a = t 2 / 3 dust 
evolution for the scale factor leaves the standard cosmo- 
logical tests unaffected. This is just as one expects for 
the very variations indicated by the observations of Jul . 




log(log t) 

FIG. 1. Numerical solution to the equations in the 
dust-dominated epoch. %p is plotted against log(logt), with 
initial conditions ip = 0,1,2,2.5. The numerical solution 
clearly approaches the asymptotic solution in the expected 
manner. The time is plotted in Planck units of 10 -43 s. 



B. The Radiation-dominated era 

In the radiation era we assume k = A = and take 
a = t 1 / 2 as the leading order solution to (0). We must 
now solve 



d_ 

dt 



Nexp[-2ip] 



There is a simple particular exact solution 

= iln(8JV) + iln(t) 
Consider a perturbation of this solution by f(t) 

^ = Jln(8A0 + iln(t) + /(t) 



(27) 



(28) 



Inserted in eqn. (27) we then get 



f +Yf f 



8t 2 



(exp [-2/]-!) 



(29) 



Let us first consider the case of a large perturba- 
tion, exp (-2/) < 1. The RHS of (||) then reduces 
to — l/(8i 2 ), and through a straightforward integration 
we get 



/ = -l + cr 3 / 2 



(30) 



with C an arbitrary constant. As t increases this will 
approach — l/(4i) which has the same absolute value and 
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is opposite in sign to the derivative of the exact solution 
(^8|) . Thus for values of tp much higher than this solution 
ip is zero, ip will stay constant until the perturbation / 
becomes small and tp approaches the exact solution (28). 



To establish the stability of the exact solution we need 
to consider small perturbations around it. For small / 
we have 



•• 3 • 1 

^ + 2t^ + Ir 2 ^ 



0. 



Hence, 



/ = -{Asin[V3Mt)] +Scos[v / 31n(t)] 



Thus, we have 



iln(8JV) + iln(t) 



+i{Asin[\/31n(i)] + B cos[V31n(i)]} 



a = 



8ATt 1/2 exp[-{Asin[\/31n(t)] 



+Scos[V31n(t)]}] -» 8Nt 1/2 



(31) 



(32) 



(33) 



(34) 



FIG. 2. Numerical solution to the equations in the radia- 
tion-dominated epoch given different initial conditions. The 
particular exact solution is eventually reached in all cases. 
The time is plotted in units of the Planck time. 



as t — > oo. 

We need to check that the ip 2 term does not dominate 
as t — > co. We have 



tp ~ 1 — 7T x oscillations 

4t t A 



(35) 



Thus the ip 2 term is the same order of t as the radia- 
tion density term if we assume a ~ t 1 / 2 . Also, the mat- 
ter density term p m (1 + C ex P [ — 2"0]) ~ Pm exp [— 2ip] ~ 
a -3 exp [— 2t/>] ~ i -3 / 2 x i -1 / 2 ~ i~ 2 is the same order of 
time variation as the radiation-density term because of 
the variation in a. The assumption a — t 1 / 2 is still good 
asymptotically but there is an algebraic constraint from 
the Friedmann eqn. (|^) 

Evaluating the terms in (^]), we have 



1 

It 2 



lirG 



M r 



S 



8m 1 / 2 



r 

t 2 



32i 2 



(36) 



where p m = Ma~ 3 , p 7 exp [-2^] = Tor 4 , iV = 2MCw~ 1 
and we have ( ~ 0.02% — 0.1% and probably ui <~ 1. So, 
to 0(t~ 2 ), we have the algebraic constraint 



8ttG 3w 
"3" ^32 



r] 



iV, are shown in fig.(|2j). The particular solution ( f28[ ) is 
clearly an attractor. It is also seen that if the system 
starts off with values higher than l/2ln(8N), ip will stay 
constant until it reaches the value of the solution, as pre- 
dicted above. In cosmological models containing matter 
and radiation with densities given by those observed in 
our universe this is the case, as seen in the computations 
shown in ref. jy. Hence, during the radiation era a re- 
mains approximately constant until the dust era begins. 

This analysis can easily be extended to other equations 
of state. If the Friedmann equation contains a perfect 
fluid with equation of state p = (7 — l)p with 7 ^ 0, 1, 2 
then there is a late time solution of (||) and (0) of the 
form 



t Si 



iV 7 2 



1)(2-7) J 



7-1 
7 



In(t) 



(37) 
(38) 



which reduces to (EST) when 7 = 4/3. This solution only 
exists for fluids with 1 < 7 < 2. 



C. The Curvature-dominated era 



This generalises the familiar general relativity (lu = 0) 
radiation universe case where we have T = 3/32nG. 

Again, the asymptotic behaviour in eqns. fl33|)-(34), 
and the approach to the exact solution (28). can be con- 
firmed by numerical solutions to eqns.(^J) - in the case 
of radiation domination. The results from runs with ini- 
tial values for ip = —8,0,8, ip = and same value for 



In our earlier study jlj] we showed that the evolution 
of a stops when the universe becomes dominated by the 
cosmological constant. This behaviour also occurs when 
an open universe becomes dominated by negative spatial 
curvature. In a curvature-dominated era we assume that 
(|) has the Milne universe solution with 
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a = t. (39) 
We must now solve eq. (^) again. It has the form 

(^t^j = ATexp[-2V>]. (40) 
We seek a solution of the form 



d_ 

dt 



4> = 2 + /(*) 

Hence, for small / 

/ + ?/ + f/ = » 

Solutions exist with / oc t n and 



-1± V1-2AT 



(41) 



(42) 



(43) 



Since N > we see that the real part of n is always 
decaying and so 



tp — ► const 
as t — > oo. Thus, as t — > oo we have 



(44) 



(45) 



where «oo and A are constants. 

Again we need to check that the tp term does not 
come to dominate. We have tp ~ £ 2 ( n_1 ) as t — > oo and 
this always falls faster than ka~ 2 oc £~ 2 since n < 0, so 
our approximation is always good. Thus we have shown 
that in open Friedmann universes a rapidly approaches 
a constant value after the universe becomes curvature 
dominated. The rate of approach is controlled by the 
matter density through the constant N in eq. (|45|) . 

This behaviour is again confirmed by numerical solu- 
tion. Fig. (||) shows how alpha changes through the dust- 
epoch and how the change comes to an end as curvature 
takes over the expansion. 



D. The Lambda-dominated era 

We can prove what was displayed in the numerical re- 
sults of [EJ, and again in fig.{Q) for the A-dominated era 
when the value of A matches that inferred from recent 
high redshift supernova observations |2^]. At late times 
we assume the scale factor to take the form 



a = exp[Xt] 
where A = and so eqn.(^|) becomes 

^(# 3At ) =Aexp[-2^] 
Linearising in tp, we have 



(46) 



(47) 




log a 

FIG. 3. The top plot shows evolution of a from radiation 
domination through matter domination and into curvature 
domination where the change in a comes to an end. The lower 
plot shows radiation (dotted), matter (solid) and curvature 
(dashed) densities as fractions of the total energy density 



Hence, 



i> + 3\Tp = Nexp[-3Xt\. 



Nt 



(48) 



ip = ip + A exp[-3A£] - — exp[-3A£] -> ip (49) 

3A 

as t — > oo, where A,ipa are arbitrary constants. Thus 
a approaches a constant with double-exponential rapid- 
ity during a A-dominated phase of the universe. The 
dominant term controlling the late-time approach to the 
constant solution is proportional to the matter density 
via the constant N. 



E. Inflationary Universes 

The behaviour found for lambda-dominated universes 
enables us to understand what would transpire during 
a period of de Sitter inflation during the early stages of 
a varying-a cosmology. It is straightforward to extend 
these conclusions to any cosmology undergoing power- 
law inflation. Suppose the varying-a Friedmann model 
contains a perfect fluid with p — (7 — l)p and < 7 < 
2/3. The expansion scale factor will increase with a(t) oc 
£ 2 / 37 , while tp will be governed, to leading order by 



(tpt 



2/7^ 







Hence, for large expansion 

tp = tp + Dr {2 -t )h -> tp Q 



(50) 



(51) 



G 



parametric form.) we have 




log a 

FIG. 4. The top figure shows numerical evolution of a 
from radiation domination through matter domination and 
into lambda domination where the change in a comes to an 
end. The lower plot shows radiation (dotted), matter (solid) 
and lambda (dashed) densities as fractions of the total energy 
density 



and so ip and a approach a constant with power-law 
(exponential) rapidity during any period of power-law (de 
Sitter) inflation. If we evaluate the kinetic term 0(ip 2 ) in 
the Friedmann equation and the terms 0(N exp[— 2ip]) in 
the ip conservation equation, we see that the assumption 
of a(t) oc i 2 / 37 is an increasingly good approximation as 
inflation proceeds. Similar behaviour would be displayed 
by a quintessence field which violated the strong-energy 
condition and came to dominate the expansion of the 
universe at late times. It would turn off the time variation 
of the fine structure constant in the same manner as the 
curvature of lambda terms discussed above. Note that 
the -0 field itself is not a possible source of inflationary 
behaviour in these models. We are assuming that the 
inflation is contributed, as usual, by some other scalar 
matter field with a self-interaction potential. However, 
if this field was charged then these conclusions could be 
altered as the coupling of the inflationary scalar field to 
the ip field would be more complicated. 



F. The Very Early Universe (t -> 0) 

As t — > we expect (just as in Brans-Dicke theory) 
to encounter a situation where the kinetic energy of ip 
dominates the evolution of a(t). This is equivalent to the 
solution approaching a vacuum solution of (^j)-(f^) with 
Pm — Pr — 0, as t — > 0. In the flat case with A = (the 
k ^ and A ^ cases can be solved straightforwardly 
and the models with p r ^ can also be solved exactly in 



AttGuj ;, 

o 



ip + 3Hip = 

Thus the exact vacuum solution is 

1 



ip = ipo + 
a = t 1 / 3 



V12ttGw 



n(t) 



(52) 
(53) 

(54) 
(55) 



During this phase the fine structure constant increases as 
a power-law of the comoving proper time as t increases: 



a 



exp[2?/>] oc t ^ 



(56) 



Note that the matter and radiation density terms fall 
off slower than ip 2 oc t~ 2 as t — > and exp[— 2ip] oc 
£-i/(V3wGw) _ They will eventually dominate the evolu- 
tion at some later time and the vacuum approximation 
will break down. As in Brans-Dicke cosmology [g6| we 
expect the general solutions of the cosmological equa- 
tions to approach this vacuum solution as t — > and to 
approach the other late-time asymptotes discussed above 
as t — > oo. 



IV. DISCUSSION 

The overall pattern of cosmological evolution is clear 
from the results of the last section even though it is not 
possible to solve the Friedmann equation exactly in most 
cases. There are five distinct phases: 

• a. Near the initial singularity the kinetic part of 
scalar field ip will dominate the expansion and the 
universe behaves like a general relativistic Fried- 
mann universe containing a massless scalar or stiff 
perfect fluid field, with a = i 1 / 3 . During this 'vac- 
uum phase', the fine structure constant increases 
as a power law in time. 

b. As the universe ages the radiation density will 
eventually become larger than the kinetic energy of 
the ip field. In this radiation dominated epoch, the 
fine structure constant will approach a specific so- 
lution, a oc i 1 / 2 asymptotically. In reality however, 
if the initial value of a is much larger than the spe- 
cific solution, we will have a potentially very long 
transient period of constant evolution, and the uni- 
verse may become dust dominated while a is still 
constant. 

c. After dust domination begins, a slowly ap- 
proaches an asymptotic solution, a = 2Aln(t) x 
exp[— f _1 Zz(t)], where U{t) is the logarithmic inte- 
gral function. If the universe has zero curvature and 
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no cosmological constant this will approach the late 
time solution a oc ln(i). 

d. If the universe is open then this increase will be 
brought to an end when the universe becomes dom- 
inated by spatial curvature and a will approach a 
constant. If the curvature is positive the universe 
will eventually reach an expansion maximum and 
contract so long as there are no fluids present which 
violate the strong energy condition. The behaviour 
of closed universes also offers a good approxima- 
tion to the evolution of bound spherically symmet- 
ric density inhomogeneities of large scale in a back- 
ground universes and will be discussed in a separate 
paper. 

e. If there is a positive cosmological constant, the 
change in a will be halted when the cosmological 
constant starts to accelerate the universe. If any 
other quintessential perfect fluid with equation of 
state satisfying p < — p/3 is present in the universe 
then it would also ultimately halt the change in a 
when it began to dominate the expansion of the 
universe. 

To obtain a more holistic picture of the evolution it is 
useful to string these different parts together. To a good 
approximation we know that in the vacuum phase from 
the Planck time t p until t v we have 

a oc t^;a oc F 4 ; A = 1 (57) 
v SttGuj 

In the radiation era we have a constant until the growth 
kicks in at a time t gr0 wth- The fine structure constant 
then increases as 

oocqoc t 1/2 (58) 

until t eq when the radiation era end and dust takes over. 
However, in universes like our own, this growth era is 
never reached. Then, in the dust era, 

a oc hit (59) 

until the curvature or lambda eras begin at t c or t\, af- 
ter which a remains constant until the present , to. So, 
matching these phases of evolution together we can ex- 
press a(to) in terms of a(t p ) : 

When the universe is open with A = : 




where we have used the fact that our log formula to 
press ages in Planck time units. 

When the universe is flat with A > : 




and t c has been replaced by t\. 

For the radiation era we consider two extreme cases. 
We look at a constant a scenario with t gr0 wth = t eq and 
a scenario where it grows throughout the radiation era, 

tgrowth — • 

Typically, t c /t p ~ t A /t p ~ 10 59 and t eq /t p ~ 10 53 , so 
in both cases for constant a evolution in the radiation 
epoch we get 

a(t ) = a(t p ) A ~ l.Ua(t p ) ^ (62) 

We approximate the value for t v *~ t p ~ 1, so for contin- 
uous growth through radiation epoch we get 

<*.>-«« (lo-r (|) ~ lo'Vfto (±y 

(63) 

Hence there are very different possibilities for the change 
in a depending on the evolution in the radiation era. 

We have proved this sequence of phases by an exhaus- 
tive numerical and analytical study. The ensuing sce- 
nario finds two interesting applications, with which we 
conclude. 

In (!]] we found that our theory could fit simultane- 
ously the varying a results reported in and the 
evidence for an accelerating universe presented in p7| . 
We noted the curious fact that there is a coincidence 
between the redshift at which the universe starts accel- 
erating and the redshift where variations in a have been 
observed but below which a must stabilise to be in ac- 
cord with geochemical evidence [p8|p9| . This may be 
explained dynamically in our theory by the fact that the 
onset of lambda domination suppresses variations in a. 
Therefore a remains almost constant in the radiation era, 
undergoes small logarithmic time-increase in the matter 
era, but approaches a constant value when the universe 
starts accelerating because of the presence of a positive 
cosmological constant. Hence, we comply with geolog- 
ical, nucleosynthesis, and microwave background radia- 
tion constraints on time-variations in a, while fitting si- 
multaneously the observed accelerating universe and the 
recent high-redshift evidence for small a variations in 
quasar spectra. 

We have also noted that within this theory the usual 
anthropic arguments for a lambda free universe may be 
reversed Q. Usually, the anthropic principle is used to 
justify the near flatness and A w nature of our uni- 
verse since large curvature and lambda prevents the for- 
mation of galaxies and stars from small perturbations. 
We have shown that it might be anthropically disadvan- 
tageous for a universe to lie too close to flatness or for 
the cosmological constant to lie too close to zero. This 
constraint occurs because "constants" change through- 
out the dust-dominated period when the curvature and 
lambda do not influence the expansion of the universe. 
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. The onset of a period of lambda or curvature dom- 
ination has the property of dynamically stabilising the 
constants, thereby creating favourable conditions for the 
emergence of structures. If the universe were exactly flat 
and lambda were exactly zero then a would continue to 
grow to a value that appears to make living complexity 
impossible |3G]. 
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